import numpy as np
from astropy import units as u

a = 45.0 * u.deg
b = np.arctan(np.tan(a))
print(b)
